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I. INTRODUCTION 

The production of top quark pairs is an interesting process for understanding QCD dy- 
namics. The very short life-time of a top quark and its large mass enable accurate theoretical 
predictions for this process for all energies, including close to ti threshold. One can then 
use the ti production to study properties of the top quark, such as the value of the top 
quark mass, branching fraction T(t — > Wb)/T tot etc. In addition, the short life-time and 
large mass of the top quark have important consequences for top quark polarization since 
non-perturbative fluctuations of chromomagnetic field are too weak to change the direction 
of the top quark spin. Hence, in the absence of hard gluon radiation, top quark polarization 
is conserved; information about it can be obtained from properties of the top quark decay 
products. This information can be used to study the Lorentz structure of interaction vertices 
involved in top quark production and decay. Also, one can use top quark spin correlations to 
distinguish different production mechanisms, e.g. qq — > ti from gg — ► ti, at least in certain 
kinematic limits. This rich physics is fully reflected in an experimental top quark program 
at the Tevatron, which is the primary source of our knowledge about top quark properties. 

n □ m 

The summary of recent experimental results can be found in Refs. [II, 12J, 13J . 

On the other hand, as was repeatedly emphasized in the past, top quark physics at the 
LHC is not just a rescaled version of the top quark physics at the Tevatron. Indeed, the 
high energy and luminosity of the LHC will lead to a dramatic increase in the number 
of observed ti events. As the result, interpretations of cross-section measurements will be 



subject to 



theoretical uncertainties rather than statistical errors. This fact motivated recent 



analysis and reviews where the quality of current understanding of tt 

production was thoroughly assessed. 

We point out in this regard that pioneering studies of heavy quark production cross- 
sections at next-to-leading order (NLO) in QCD were performed almost twenty years ago 
js, 9] 1 . These computations were furt 
the produced top quarks in Refs. [U 



rer extended to cover various kinematic distributions of 
121 ] . In addition, threshold resummation was applied to 



top quark production at the Tevatron and the LHC 13], although its practical relevance for tt 



production at both colliders is still an issue of debate. Soft gluon corrections and Coulomb 



It is interesting to point out that analytic results for NLO QCD corrections to top production in gg, qq 
and qg channels were obtained very recently in Ref. [To| . 



eluon bound state corrections were combined in [l4| to provide an accurate description 

n 

of tt threshold region. Electroweak corrections to a t i were studied in [15(. There is an 
ongoing effort to compute next-to-next-to-leading order QCD corrections to top quark pair 



production in hadron collisions 161 ]. 



In all studies of top production described so far top quarks were treated as stable particles 
and summation over their spin degrees of freedom was performed. This approach introduces 
two problems. First, spin degrees of freedom of top quarks influence kinematics of their decay 
products and in this way lead to observable consequences. Second, there are QCD radiative 
corrections related to the decay, rather than production, stages of the time evolution of ti 
system. Given the fact that the tt production cross-section at the LHC is very large, a 
degree of realism in its description is clearly warranted. To achieve it, we require a NLO 
QCD prediction for ti production that is valid at the level of observable particles, such 
as leptons, quarks and gluons originating either from production of top quarks or in their 
decays. In principle, this seems to necessitate a next-to-leading computation for 2 — > 4 
processes such as pp — * l + vi~i)hb or pp — > ud£~vbb, which is a formidable task at present. 

Fortunately, the problem can be simplified by studying double resonance contributions, 
but accounting for spin degrees of freedom exactly through all stages of the top quark decay 
chain. Indeed, when top quarks are treated as truly unstable particles, all QCD corrections 
to a relevant 2^4 process 2 can be decomposed into factorizable and non-factorizable jl3] ■ 
Non-factorizable corrections imply a cross-talk between production and decays of top quarks. 
In the limit Y t /m t — > 0, these corrections must vanish since quantum interference should 
not occur if events are separated by macroscopic distances. The precise way in which such 
non-factorizable corrections vanish was described in Refs. \v\ , [la . . In what follows we 
work in the on-shell approximation for top quarks, ignoring non-factorizable corrections. 

Once non-factorizable corrections are neglected, a full description of tt production and 
decay including all the spin correlations is achieved by computing NLO QCD corrections to 



both production and decay of a polarized ti pair. In an impressive series of papers |20l . l2lL 122 ] 
Bernreuther, Brandenburg, Si and Uwer computed the spin density matrix for the production 
of a ti pair in hadron collisions through NLO QCD. Corrections to decays of polarized top 



quarks were obtained in Refs. 23j, |24j. Putting all these bits together, Bernreuther et al. 



2 We mean here the hnal state W + W bb. 
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studied a number of kinematic distributions which can be used at hadron colliders to probe 



top quark spin correlations 



25 



While the importance of accurate predictions of top quark spin correlations was strongly 
emphasized in Refs. 25|, |26|], to the best of our knowledge there is no publicly available 
numerical program that satisfies the following requirements: 

• it containts NLO QCD corrections to top quark production and decay; 

• it includes all spin correlations for the processes pp — ■> £ + v£~vbb and pp — > udi^ubb; 

• it allows arbitrary cuts on particles in the final state. 



While the implementation of heavy flavor production in MC@NLO [27] and POWHEG [28] 
comes close to these requirements (see e.g. 29]), those programs do not fully include all 
spin correlations through NLO QCD. It seems to us that developing a numerical program 
with capabilities listed above is useful since it will contribute to more realistic description of 
ti production at hadron colliders. Computing NLO QCD corrections to ti production and 
decay and implementing them into a flexible numerical program is the primary goal of this 
paper. 

Another motivation for undertaking the study of ti production is more theoretical - 
we would like to explore how unitarity-based methods for one-loop computations work in 
a relatively simple but fully realistic setting, when massive particles are involved. Initial 
studies of ti pro duction in the context of generalized D-dimensional unitarity 30| were 
performed in 3jJ- Those results were extended by us with an eye on applying generalized 
D- dimensional unitarity to ti production in association with jets 3 . The simplest application 
- the case of polarized ti production - is described in this paper. 

The rest of the paper is organized as follows. In Section [TT] some theoretical aspects of 
the computation are discussed. In Section II III we present a number of observables relevant 
for ti studies at the Tevatron and the LHC. We conclude in Section HVl More details on the 
calculation of virtual and real corrections are described in Appendices [A] and UU respectively. 



3 We point out that NLO QCD corrections to tt+jet production at the LHC and the Tevatron were reported 
in 0,1331. 
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II. THEORETICAL FRAMEWORK 

In this Section the theoretical framework relevant for the computation is briefly described. 
As already pointed out in the Introduction, we study the production of top quarks and allow 
for their semileptonic decays. We include all the spin correlations but keep top quarks on 
their mass shells. We do not consider hadronic decays of top quarks in this paper. 

We first discuss an efficient way to incorporate top quark decays into the computation 4 . 
The tree amplitude for the process ij — » it — > (M~u)(b£ + u) can be written as 

A«~ = U{t - bfu) 2 i( * 2 +m O A(ij - ft) ( 2 l{ -^ 2 +mt \ i(t - 5TP 
\ pf — ml + rni t L t I \ p| — mf + im t i t 

where A(ij — >i£) and A(t — > 6 £ ± ( z/') are sub- amplitudes for production and decay pro- 
cesses, top quark propagators are factored out and summation over spinor indices is implicit. 
Note that in Eq.([T]) the top quarks are off the mass-shell. To compute the cross-section, 
we need to square the amplitude y4 tree and integrate it over the phase-space for final state 
particles. To simplify this procedure, we approximate the squared amplitude by taking the 
limit T t /m t — > 0. In this approximation, propagators that appear in Eq.([T]) yield 

277 8(pl-mt). (2) 




{P 2 t -m 2 t ) 2 + m\Y\ 



2m t T 



After factorizing the phase-space integral into the phase-space for a ti pair and the decay 
phase-spaces for t and i, the delta- function in Eq.p]) forces top quarks on their mass-shells. 
This separates the production stage from the decay stage and allows us to implement top 
quark decays by choosing the on-shell spinors 

U( Pt ) = A(t - b£ + v) (3) 

V(p,)='tpplA(t^M-, ) , (4) 

V2m t r t 

to describe polarization states of top quarks. Eq.(pQ) can now be rewritten as 

A trec = U(pt) A(ij - it) V( PI ) + O (J^j . (5) 



We are grateful to Keith Ellis for emphasizing this point to us. 
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The usefulness of Eqs.([3])-(j5]) is twofold. First, since U and V are on-shell Dirac spinors, 
the amplitude A tree can be computed in a conventional way; this means that the inclusion 
of t and i decays does not increase the complexity of scattering amplitudes that need to 
be calculated. Second, Eq.Q actually helps in reducing computational burden. Indeed, 
the top quark spinors depend on polarization states and momenta of leptons, neutrinos and 
6-quarks which we treat as massless. For top quark decays, helicity states of all massless 
particles are fixed, due to the V — A structure of flavor-changing interaction vertices in the 
Standard Model. This implies that the spinors U and V have uniquely defined polarizations 
and no helicity sums related to t and t are involved in the cross-section computation. In the 
evaluation of the decay amplitudes we keep the IV-boson on-shell and treat the 6-quark as 
massless particle. 

Having discussed how to deal with t and i decays efficiently, we focus on other compu- 
tational details. The on-shell approximation for top quarks splits all the QCD effects into 
corrections to the production and corrections to the decay. These corrections do not inter- 
fere since production and decay stages are separated by large space-time intervals. As the 
result, ti production with NLO QCD effects and t(i) decay with NLO QCD effects can be 
treated as independent processes. 

We begin with corrections to it production process. In this case, as we just discussed, all 
one has to do to incorporate top quark decays, is to use in Eq. (j5]) the on-shell spinors for t 
and i, shown in Eqs.([3])-(j4]). This prescription is valid both at leading and at next-to- leading 
order. We organize the calculation in terms of gauge-invariant color ordered sub-amplitudes. 
Each such amplitude multiplies a specific color factor and has a fixed ordering of external 
particles. We use Berends-Giele [34J recurrence relations to evaluate those sub-amplitudes for 
each contributing helicity configuration at leading order. There are two partonic processes 
that need to be considered for LO, and virtual component of NLO computations - gg — > ti 
and qq — > ti. Their color decomposition is given in Appendix [A] 

For the computation of virtual corrections to ti production, we employ the method of 
generalized Z)-dimensional unitarity 30]. The basic idea of this method is to consider 
unitarity cuts of scattering amplitudes in D-dimensional space-times where D > 4 is integer. 
It turns out that dimensionally-regularized one-loop amplitudes can be fully reconstructed 
from these unitarity cuts provided that tree-level on-shell scattering amplitudes for complex 
external momenta, evaluated in six- and eight-dimensional space-times, are available. For a 
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detailed description of the method see Ref. [301 ] . 

We turn to the discussion of real emission corrections to ti production. Real emission 
processes develop singularities when a massless final state particle is unresolved, i.e. it is soft 
and/or collinear to another particle. When integrated over the unresolved phase-space, those 
singularities produce divergencies that cancel against divergencies in virtual corrections and 
unrenormalized parton distribution functions, yielding finite cross-sections. In order to treat 
real emission singularities in a numerically stable way we employ an extension of Catani- 



Seymour dipole subtraction scheme [35| to massive particles |36j| . Recall that the basic idea 
of the subtraction procedure is to construct an approximation to matrix elements squared 
of real emission processes that, on one hand, has the same singular limits as real emission 
matrix elements and, on the other hand, can be analytically integrated over the phase-space 
of unresolved particles. Because top quarks are heavy, they never appear as unresolved 
particles and their spin degrees of freedom are not essential for the construction of the 
subtraction terms. Hence, we can also use the U, V spinors to compute (subtracted) real 
emission corrections and in this way account for decays of top quarks exactly. The color 
decomposition of amplitudes needed for the real corrections to the ti production process 
is given in Appendix |A] Furthermore, a list of all dipole subtraction terms is given in 
Appendix [B] 

We now turn to the discussion of NLO QCD corrections to the top quark decay. Cor- 
rections to decay rate and lepton kinematic distributions in decays of polarized top quarks 
are well-known {23], • However, since we want to keep our computation completely dif- 
ferential, to allow for arbitrary cuts on the final state particles, we need to go beyond the 
computation of Ref.j23|]. A detailed discussion of how such a computation should be set up 
was given recently for the case of single top production in Ref. {38] and we closely follow that 
reference in our implementation of virtual and real QCD corrections to top quark decays. 

We also point out that we use and Tf LO to construct top quark spinors in Eqs. (G])-(J4|) 
in LO and NLO computations, respectively. Indeed if no restrictions on final state particles 
are applied, the cross-section computed in the on-shell approximation is given by the product 
of the ti production cross-section multiplied by t and i decay branching fractions. Since this 
statement holds true to all orders in perturbation theory, the width in Eqs.([3~])-(jll) should 
also be computed in series of a s , for consistency. Hence, the complete NLO cross-section 
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can schematically be written as 



Ayo Aro HP- dF° 

da NLO = da ^4^- + da 1 d * d * 



(r°) 2 (r°; 



2 



+ dff « (M + 52^T) _ ia ^Il^m (6) 

+ V (r») 2 + (r») 2 J r» (r«) 2 w 

where cr°, T and a 1 , T 1 denote leading and next-to-leading order contributions to the 
production and decay processes, respectively. The last term in Eq.([6]) arises from the 
correction to the top quark width in Eqs.([3])-(j4]). 

To conclude this Section, we mention checks that we applied to our computation to 
ensure its correctness. Within the unitarity method, one computes scattering amplitudes 
for particular polarization states of external particles. As the result, gauge invariance of 
all one-loop amplitudes involving external gluons can be easily checked by substituting 
polarization vector of a gluon by its momentum. One can also check the threshold s —>■ 4mf 
behavior of one-loop amplitudes, relevant for ti production. Indeed, close to ti threshold, 
scattering amplitudes are dominated by the Coulomb singularity and we checked that our 
implementation of the virtual corrections reproduces the Coulomb singularity correctly. The 
cancellation of infra-red and collinear divergencies between virtual and real corrections was 
checked numerically. Furthermore, the implementation of both virtual and real corrections 
was checked by comparing results for physical observables produced by our code against 
other programs that compute the ti production cross-section without top decays. We have 



extensively used MadGraph [39( and MCFM [40| for these comparisons to check LO and 
NLO results, respectively. In order to ensure the correctness of our implementation of top 
quark decays we checked that the NLO QCD corrections reproduce known results for the 
top width if no restrictions on final state particles are applied. Finally, we have checked that 
our code reproduces NLO QCD corrections to some observables sensitive to top quark spin 
correlations, discussed in Ref.[26j. 



III. RESULTS 



We have implemented NLO QCD corrections to the production and semileptonic decays 
of a top quark pair in a FORTRAN program. The goal of this Section is to illustrate 
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some of its potential applications. Before discussing those examples, we describe the input 
parameters that are used in all numerical results presented in this paper. 

We consider production of ti pairs both at the Tevatron and the LHC. The Tevatron 
center of mass energy is yfs = 1.96 TeV and the LHC center of mass energy is y/s = 10 TeV. 
We use the pole mass of the top quark, mt = 172 GeV. The strong coupling constant 
as is renormalized in the MS scheme except for top quark loop contributions to gluon 
self-energy diagrams that are subtracted at zero gluon momentum. We employ CTEQ6L1 
parton distribution functions for leading order and CTEQ6.1M for next-to- leading order 
computations [4l|, |42J]; those PDF sets correspond to a s (rriz) = 0.130 and ot s {mz) = 0.118 
respectively. The mass of the W is taken to be mw = 80.419 GeV. W couplings to fermions 
are obtained from the Fermi constant Gf = 1.16639 • 10~ 5 GeV -2 . The leading order width 
of the top quark is = 1.47 GeV and Tf LO = 1.31 GeV at next-to-leading order 5 . We use 
|Vtb| = 1 and employ on-shell approximation for the W boson, produced in t — > Wb decays. 

We set the width of the W boson to 2.14 GeV. We consider this to be an experimentally 
measured width of the W boson, not the result of computation in leading order of QCD 
perturbation theory. This implies that our leading order prediction for pp — > ti — > l^vi'vbb 
includes the NLO QCD branching ratio for W — > ev. While simple, this effect is definitely 
not negligible numerically - since 0(a s ) corrections increase the W width by approximately 
2.5 — 3 percent, the use of LO or NLO branching fractions for W bosons changes the cross- 
section for pp — > ti — > i^vt^vbb by five to six percent. 



To define jets, we use the /^-clustering algorithm |43J] with R = 0.4. We require that 
two 6-jets are present in the event. We define jet flavor through its "bottomness" quantum 
number so that, for example, a jet that contains b and b is not a 6-jet. We point out that this 
definition is infra-red safe through NLO QCD approximation for ti production whereas this 
is not true in general, for massless quarks 44j. For 6-jets we require a minimal transverse 
momentum of 20 GeV. Charged leptons must be produced with the transverse momentum 
larger than 20 GeV and the missing energy in the event should exceed 40 GeV. We also 
set renormalization and factorization scales to the value of the top quark mass. Since scale 
dependence of various observables for ti production was studied in detail in the existing 



5 This result for rf LO is obtained by choosing the renormalization scale for the strong coupling constant 
to be the top quark mass. 
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FIG. 1: Various kinematic distributions for leptonic final states in tt production at the Tevatron. 
We show transverse momentum (a) and rapidity (b) distributions of the positively charged lepton 
as well as the distribution in the invariant mass of a lepton and a 6-jet (c). The distribution in the 
(particularly defined, see text) opening angle of the leptons is shown in (d). Each distribution in 
panel (d) is normalized to the corresponding total cross-section. All cuts described at the beginning 
of Section III are applied. 

literature, we do not present such studies in this paper. 

To set the scale for the magnitude of next-to-leading QCD effects in tt production for our 
choices of input parameters, we quote results for cross- sect ions at leading and next-to-leading 
order for the Tevatron (pp — > tt — > bl + vbl~v) 

C"NLO 



olo = 34.63 fb, cr NLO = 36.47 fb, K TEV 



1.05, 



and the LHC (pp tt -> bl + vbl~v) 



a LO = 1484fb, a NLO = 2097fb, K LRC 



C"LO 



0"NLO 
0"LO 



1.41. 



(7) 



(8) 



To obtain those numbers, we set the renormalization and factorization scales to m t and 
apply all the cuts listed in the beginning of this Section. 
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FIG. 2: Various kinematic distributions for leptonic final states in tt production at the LHC. We 
show transverse momentum (a) and rapidity (b) distributions of the positively charged lepton as 
well as the distribution in the invariant mass of a lepton and a 6-jet (c). The distribution in the 
(specially defined, see text) opening angle of the leptons is shown in (d). Each distribution in panel 
(d) is normalized to the corresponding total cross-section. All cuts described at the beginning of 
Section III are applied. 

We are now in position to illustrate capabilities of our numerical program by presenting 
a number of tt-related kinematic distributions, computed through NLO in perturbative 
QCD. In Fig{T]we present results for the Tevatron. The transverse momentum and rapidity 
distributions of leptons in top decays are shown in Figs. [T^a) andQJb), respectively. The 
distribution in the invariant mass of the charged lepton and the 6-jet is given in Figfl^c). 
Figfjjd) shows the distribution in cos ipi+i-, where (fi+i- is the angle between the directions 
of flight of £ + and £~, defined in the rest frames of t and i respectively. In all cases we 
compare predictions at leading and next-to-leading order, with and without corrections to 
the decay. Corresponding results for the LHC are shown in FigfSJ 

We first consider transverse momentum and rapidity distributions of the charged lepton, 
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FIG. 3: The ratio of NLO to LO predictions for £ + transverse momentum distributions, for the 
Tevatron and the LHC. The input parameters are described at the beginning of Section III. 

shown in the upper panels of Figs JT1l2l Both of these distributions are standard but effects of 
QCD corrections to top decays and top spin correlations are never included in their compu- 
tation. For both of these observables NLO QCD effects are important. For example, as we 
show in FigEJ the shape of the transverse momentum distribution changes significantly at 
the Tevatron while at the LHC similar change is smaller, but non-negligible. For the charged 
lepton rapidity distribution, NLO QCD corrections to decays are important, especially at 
central rapidities at the Tevatron. 

Another distribution that we show in Figs JTO is the distribution in the invariant mass of 
a positively charged lepton i + and a 6-jet, M^+&. This observable is interesting for a number 
of reasons. First, a (kinematically) similar observable - an invariant mass of £ + and J/ip 
originating from B meson decay - was discussed in connection with the measurement of the 
top quark mass in Ref.{45]. Similar to J ftp £ + invariant mass, the invariant mass of the 
6-jet and i + has a clear kinematic boundary at leading order, see Figs JTf2l This boundary 
max(M £ 2 + J = mj — m 2 ^ is the consequence of the on-shell conditions for the top quark and 
the W boson. Those conditions lead to a lower bound on the allowed neutrino energy in the 
top quark rest frame E u and to the upper bound on Mf +b = m 2 — 2m t E u . The existence 
of the kinematic boundary makes this observable potentially interesting for the top quark 
mass determination. In addition, Mi+ b belongs to a general class of observables that may be 
used to study spins of Beyond the Standard Model particles at the LHC. The effect of NLO 
QCD corrections on the M^+f, distribution is interesting. For example, at the Tevatron there 
is a significant cancellation between corrections to ti production and corrections to t and i 
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decays for this observable. This effect also exists at the LHC although it is less pronounced. 

It is interesting that once the NLO QCD corrections are included, events start to appear 
beyond the leading order kinematic boundary and both, QCD corrections to the production 
and QCD corrections to the decay contribute. The origin of these events is peculiar since 
they appear as the consequence of the fact that the momentum of the fe-jet rather than 
the momentum of the 6-quark is employed in the calculation of M^. Then, for example, 
a gluon emitted in the decay of i may combine with the 6-quark from top decay to form 
a 6-jet. The energy of such a jet is not restricted by on-shell conditions applied to t and 
W. As the result, the invariant mass of the 6-jet and £ + may exceed the LO kinematic 
boundary. A similar mechanism is in effect when gluon emission occurs as part of the top 
production process. Note, however, that if a gluon is emitted in the decay of a top quark 
t — > hi^v + g, then the invariant mass of a 6-jet and £ + can not exceed the leading-order 
kinematic boundary. 

Finally, we discuss lepton angular correlations; such observables are particularly sensitive 
to correct implementation of spin correlations in ti production and decay. Indeed, angular 
correlations tell us whether leptons prefer to be produced with parallel or anti-parallel 
momenta and, as it turns out, the answer to this question differs for the Tevatron and the 
LHC [7]. To understand the difference recall the two mechanisms that dominate the ti 
pair production at those colliders. At the Tevatron, top pairs are mostly produced in an 
annihilation of a qq pair - qq — > g* — > ti - which forces a ti pair to have angular momentum 
J = 1. Since, in addition, most of the time the annihilation occurs close to ti threshold, it 
is most probable that the ti pair is produced in an .S-wave with spins of t and i parallel, 
to create a J = 1 state. Since e + likes to follow the direction of t spin while e~ prefers the 
direction opposite to i spin, flight directions of e + and e~ are anti-correlated. 

At the LHC the situation is different since gg — > ti annihilation becomes the dominant 
production channel. Close to ti threshold the largest contribution comes from an S"-wave 
J = color-octet annihilation 14j . This suggests that e + and e~ prefer to be produced with 
parallel momenta at the LHC, at least to an extent that threshold production is important 
there. 

In Figs. [U(d)j2](d) we show the opening angle distribution between two charged leptons 
at the LHC and the Tevatron. Note that the flight directions of i + and l~ are defined, 
respectively, in t and i rest frames, when computing this observable. Azimuthal correlations 
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FIG. 4: The distribution in the opening angle of the two leptons in the laboratory frame, normal- 
ized to the corresponding total cross-section. All cuts described at the beginning of Section III are 
applied. 

are substantial at LO and remain fairly pronounced even after NLO QCD corrections are 
included. It is interesting that NLO QCD effects are more important at the Tevatron where 
the shape of the distribution changes. At the LHC the NLO QCD effects do not change 
the shape of the distribution at all, so leading order predictions do a good job in that case. 
QCD corrections to top decays do not play an important role for this observable for both, 
the Tevatron and the LHC. We note that a similar distribution was computed through NLO 
QCD in Refs. (25I 26]. but no jet cuts were applied. However, by removing all the cuts 
in our computation and by adjusting the input parameters, we reproduce results for this 



distribution reported in Refs. 



25 



26) 



One of the drawbacks of using cos tpi+i- to study spin correlations is the requirement that 
t and i rest frames are reconstructed and, for semileptonic top decays, it is not possible to 
do this unambiguously. In this respect, we would like to point out that ti spin correlations 
can also be studied using a much simpler observable - the opening angle ipe+e- of the two 
leptons in the laboratory frame. The corresponding distributions computed at leading and 
next-to-leading order are shown in FigJH for both the Tevatron and the LHC. The results 
are similar to distributions in <f£+£- shown in Figs. []Jd)j2^d) but the effect of the radiative 
corrections on the shape of ipe+e- distribution at the LHC is stronger. 

To conclude this Section, we point out that results reported here are not supposed to 
provide exhaustive phenomenological studies of ti production; rather, they should illustrate 
capabilities of our numerical implementation of NLO QCD corrections to ti production and 
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decays. We have chosen to present a number of distributions which we find interesting but, as 
with any NLO QCD computation, essentially any infra-red safe observable can be calculated 
and analyzed. In this regard, a numerical program that allows to impose arbitrary cuts and 
employ arbitrary jet algorithms, includes all the spin correlations in tt production and in t 
and i decays and incorporates NLO QCD corrections to (semileptonic) decays of top quarks, 
is very attractive since it can be used for realistic description of tt pair production in hadron 
collisions. 



IV. CONCLUSIONS 

In this paper the computation of next-to-leading order QCD corrections to the production 
and decay of top quark pairs at the Tevatron and the LHC is presented. Our goal was to 
develop a numerical program which employes the on-shell approximation for top quarks, yet 
accounts for all the spin correlations in their production and semileptonic decays through 
NLO QCD. 

To this end, the implementation of next-to-leading order corrections to the production 
and decay of polarized top quarks is required. We performed such an implementation at a 
fully differential level. This allows us to compute an arbitrary observable defined in terms 
of lepton momenta, missing energy and jet energies and momenta, originating either in top 
production or decay. We have illustrated the capabilities of the program by presenting a 
variety of differential distributions that are sensitive to top spin correlations and NLO QCD 
corrections to top quark decays. 

An interesting aspect of the computation reported in this paper is that the method of 



generalized .D-dimensional unitarity [30|] is employed to compute virtual corrections to tt 
production cross-section; this is the first application of generalized .D-dimensional unitarity 
to a fully realistic NLO computation that involves massive external and internal particles. 
We find that, in general, the method works well so its application to other processes with 
massive particles is definitely warranted. 
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APPENDIX A 

In this Appendix we summarize various results for the color decomposition of tree- and 
one-loop amplitudes that we used in this paper. 

We begin with the leading order process. At leading order two partonic processes con- 
tribute to ti production: gg — > ti and qq — > ti. Their color decomposition is given by 

A^(gg -> ti) = gl £ (T"-»T-*)£ ^tree (1 _ ^ ^ (Al) 

<T&S 2 

A*"{qq - tt) = g 2 s - 2 <> ^ ^ ^ 

where A tree (lt,2 t ,i, j) are color-ordered tree amplitudes and g s is the strong coupling con- 
stant. The generators T a of the SU(Af c = 3) color group are normalized to Tr(T a T fe ) = 5 ab 
and satisfy the commutation relation [T a ,T b ] = —F£ b T°. 

Considering real emission contribution to top quark pair production at next-to-leading 
order, we find that four partonic processes contribute gg — > itg, qq — > itg, qg — > itq, qg — > 
ttq. The last three processes are related by crossing symmetry. The color decomposition for 
the two master processes is given by 

A *»(gg -> tig) = g 3 s ]T ( T ^ T ^T^)| A *~( lh 2t , (a 3 ) g , (a 4 ) g , (a 5 ) g ), (A3) 



A*"{qq - tig) = g 3 s [ (T^)} 4 8£ A^{l h 2 U 3„ 4„ 5 g) 



+ (r^)|^^ tree (l f ,2 t ,5 9) 3 5 ,4 g ) 
+ ^ c (T a nZSilA t ^(l,,5 g ,2 t ,3 q ,4 q ) 

+^(T as )g Sll ^ ree (l f ,2 t ,3„5 9 ,4 g )]. (A4) 
In the case of virtual amplitudes, a decomposition into color-ordered amplitudes is insuffi- 



aes 2 



cient; instead, one has to consider the so-called primitive amplitudes [461 ] . We need one-loop 
amplitudes for gg — ► ti and qq — > ti. For gg — > ti the color decomposition reads 

A™\gg -> tt") = g A s ^ [ (T^)l(F a ^F a ^) XlX2 A L ^(l h 2 t , (a 3 ) g , (a 4 ) g ) 

+(T*>T°*T x %(F a "*) xlX2 A L ^(l h (a 3 ) g ,2 t , (a A ) g ) 
+ {T X ^T^T X % A h > [l \l h (a 3 ) g , (a,) g ),2 t ) 

JVf 

+ ^(T-3 T -0| 4' [1/2l (lf,2 t >3) 3 > 4 ) 5 )] (A5) 
/=i 
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FIG. 5: Definition of primitive amplitudes with four quarks. Fig. a) defines A a which includes 
all topologies where both external fermion lines enter the loop; Figs, b) and c) define A b and A c 
where either the top quark or the light quark line enter the loop, respectively; Fig. d) A d includes 
topologies with a closed fermion loop. 



where A L '^ are left-primitive amplitudes, defined in 46[]. Nf is the number of quark flavors; 
top quarks in closed fermion loops are treated as massive particles. 



For the virtual corrections to the process qq — > tt we use the color decomposition 

1 



A ™*(qq -> tt) = g* 



N c 12 u 



(A6) 



As the second term in Eq. (IA6l) vanishes after interference with the leading order amplitude 
we do not consider it here. The first term in Eq. (1A6l) is given by 

Bi = (N c - J,)^ a (l f ,2 t ,3 5 ,4 ? )-|,^ tt (l t -,2 t ,4 ?) 3 9 -) 

+— (A\lt, 2 t , 3 q , 4 q ) + A C (U, 2 t , 3 ? -, 4,)) - ^ A d f {l h 2 t , 3 ? , 4,). (A7) 
f=i 

The ordered primitive amplitudes A a,b ' c ' d correspond to different topologies depending on 
which fermion lines enter the loop, cf. Fig. 5. 
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APPENDIX B 



Following the notation of Ref.[35[], the generic form of the subtraction terms for the 
partonic channels g(l) g(2) — > i(3) t(4) g(5) is given by twelve dipoles 

V\\ Vf, Vf, Vf, P 15 ' 2 , P 25 ' 1 , Vl 5 , V 2 35 , Vl 5 , D 2 5 , V 35A , P 45j3 . (Bl) 

It is easy to see that identical list of dipoles is required to construct the subtraction terms for 
q(l) q{2) — > i(3) t(4) g(5) partonic channel. Note, however, that the above dipoles differ for 
different partonic channels since they contain particular splitting kernels and reduced tree 
amplitudes. The processes V(l) g(2) — > t(3)t(4) q (5) develop only collinear singularities 
making a summation over all spectator particles unnecessary. We therefore choose 

£> 15 ' 2 , P 25 ' 1 (B2) 

as subtraction terms which differ from the original dipoles Eq. (IBip only by their color factor. 
In fact, since collinear singularities are local in color space, no color correlations are required 
for the dipoles in Eq. flB2l) ; the corresponding color factors are given by Cp and T R for each 
dipole in Eq. (1B2l) . respectively. 

The auxiliary contribution from the dipole subtraction terms needs to be added back. 
We choose to implement those terms by separately integrating each dipole over the unre- 
solved phase-space. Singularities of integrated dipoles cancel against singularities of virtual 
corrections and unrenormalized parton distribution functions, yielding finite result for the 
cross-section. For the implementation of both unintegrated and integrated dipoles we used 



results in Refs. 



36 



. Following Refs. 



3£ 



47 



48J we introduce a parameter a for all 



initial-state emitter dipoles which allows to cut off the dipole phase-space in non-singular 
regions. This helps to improve numerical stability of our code and provides a non-trivial 
check of the implementation of subtraction terms. 
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